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ABSTRACT 

We present new [SII] images of the HH 30 jet and counterjet observed in 2006, 2007, and 2010 
that, combined with previous data, allowed us to measure with improved accuracy the positions and 
proper motions of the jet and counterjet knots. Our results show that the motion of the knots is 
essentially ballistic, with the exception of the farthest knots, which trace the large scale "C" -shape 
bending of the jet. The observed bending of the jet can be produced by a relative motion of the HH 
30 star with respect to its surrounding environment, caused either by a possible proper motion of 
the HH 30 star, or by the entrainment of environment gas by the red lobe of the nearby L1551-IRS 
5 outflow. Alternatively, the bending can be produced by the stellar wind from a nearby CTTS, 
identified in the 2MASS catalog as J04314418+181047. The proper motion velocities of the knots 
of the counterjet show more variations than those of the jet. In particular, we identify two knots 
of the counterjet that have the same kinematic age but whose velocities differ by almost a factor of 
two. Thus, it appears from our observations that counterjet knots launched simultaneously can be 
ejected with very different velocities. We confirm that the observed wiggling of the jet and counterjet 
arises from the orbital motion of the jet source in a binary system. Precession, if present at all, is of 
secondary importance in shaping the jet. We derive an orbital period r G = 114 ± 2 yr and a mass 
function m/zj? = 0.014±0.006 M Q . For a mass of the system of m = 0.45±0.04 Mq (the value inferred 
from observations of the CO kinematics of the disk) we obtain a mass rrij = 0.31 ± 0.04 M Q for the 
jet source, a mass m c = 0.14 ± 0.03 Mq for the companion, and a binary separation of a — 18.0 ± 0.6 
AU. This binary separation coincides with the value required to account for the size of the inner hole 
observed in the disk, which has been attributed to tidal truncation in a binary system. 
Subject headings: ISM: Herbig-Haro objects — ISM: individual (HH 30) — ISM: jets and outflows — 
stars: formation — 



1. INTRODUCTION 

Th e Herbig-Haro (HH) object 30 (|Mundt fe Fried! 
1983), located in the northeastern part of the L1551 
dark cloud, lies at a distance of 140 pc (jKenvon et al.l 
1994). The HH 30 outflow is considered a prototypical 
jet/disk system. It presents a clear iet/counterjet struc- 
ture, whic h has been described by IMundt et al.1 (|1987l 
\Em> and IGraham fc Heverl IpM ). The HH 30 excit- 
ing source is an optically invisible star (iVrba et al.lll999D 
highly extinguished by an e dge-on disk ((Burrows et al.l 
Il996t iStapelfeldt et al.lll999D . which extends up to a ra- 
dius of ~ 250 AU perpendicularly to the jet, and divides 
the surround i ng re flection nebulosity into two lobes. 
iKenyon et al.l (1 19981) propose a spec tral type M0 for the 
HH 30 star, and lCotera et al.l (| 20011 ) estimate a bolomet- 
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ric luminosity of 0.2-0.9 L Q . lLopez et al.l (|1995l I1996D 
propose that a number of knots located to the northeast 
of the HH 30 object are also part of the same flow, result- 
ing in a total angular size of ~ 500" (~ 0.35 pc) for the 
whole outflow. Several studies have explor ed the spatial 
morphology along an d across the jet axi s (IMundt et ah 
1991 1 : iRav et~1dll996l) . the 3-D structure (lEsauivel et al. 
2007t lDe Colle et alJ 120101) . line ratios (IMundt et al. 
1990: iBacciotti et al] [19991: lHartigan fc Morsel I2007D . or 
radial velocities (IRaga et al.lll997f l of the HH 30 flow. 

In addition, some studies have been carried out study- 
ing the variability of the re fl ection nebulosity around HH 
30. iWatson fc Stapelfeldtl (|2007t) find variations of the 
lateral asymmetry of the nebula and counternebula, but 
find no convincing evide nce for any period. However, 
iDuran-Rojas et aLl (|2009f ) find a periodic modulation of 
the polarization of the nebula, with a period of 7.5 days, 
which is interpreted as produced by asymmetric accre- 
tion hot spots on the star, or orbiting clumps or voids in 
the disk. 

Proper motions of a few knots of t h e HH 30 flow 
were d eterm ined bv IMundt et al.l (|1990l ) , iBurrows et al.1 
(| 19961 ) . and lLopez et al.1 (|1996[ ). A thorough study of 
the proper motions of all the knots of the HH 30 jet 
and the closest knot s of the counterjet was made by 
lAnglada et aLl (|2007l ). The main result of this study is 
that the overall structure of the HH 30 jet can be well fit- 
ted by the shape of a wiggling ballistic jet, arising either 
by the orbital motion of the jet source or by precession 
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of the jet axis because of the tidal effects of a companion. 
In the first case the binary separation is expected to be 
9-18 AU, while in the case of precession the binary sepa- 
ration is < 1 AU. Given that the radius of the flared disk 
observed in the HST images is 250 AU, the conclusion is 
that this disk appears to be a circumbinary disk rather 
than a circumstellar disk, contrary to what was initially 
thought. This circumbinary disk is unlikely to have a 
r elevant role i n the jet collimation. 

iPetv et all (|2006( 1 carry out PdBI observations of the 
dust continuum and CO line emission revealing the pres- 
ence of an asymmetric molecular outflow and a disk in 
Keplerian rotation, deriving a central stellar mass o f 
0.45 ± 0.04 M Q . More recently, iGuilloteau et al.l (|2008[ ) 
carry out CO and dust continuum observations of the 
circumbinary disk that reveal an inner cavity with a ra- 
dius of 37 ± 4 AU, and lMadlener et al.l ((2012), throught 
a detailed modeling, conclude that the disk has an in- 
ner dep letion zone with a simi lar radius os 45 ± 5 AU 
radius. IGuilloteau et al.l (|2008[) explain the size of the 
inner hole as a result of tidal forces induced by a binary 
with a separation of 18 ± 2 AU, supporting the binary 
interpretation proposed by Anglada et al.l ()2007l ) . 

In this paper we present new [SII] images of the HH 30 
jet and counterjet carried out in 2006, 2007, and 2010, 
which allowed us to measure the position and proper mo- 
tion of the jet and counterjet knots, and to constrain the 
parameters of the binary system. In §2 we describe the 
observations; in §3 we describe the procedure used for 
the proper motions determination; in §4 we study the 
large-scale jet bending, and its possible origin; in §5 we 
analyze the jet wiggling and the constraints on the phys- 
ical parameters of the binary system at the core of HH 
30; finally, in §6 we discuss the results obtained and in 
§7 we give our conclusions. 

2. OBSERVATIONS 

The CCD observations used in this paper to deter- 
mine the proper motions of the HH 30 jet/counterjet 
system are listed in Table [U All the images were ob- 
tained through [SII] narrow-band filters, which included 
the AA6716, 6731 A emission lines. The details on the 
setup configuration, acquisition and treatment of the 
2.5 m Isaac Newto n Telescope (INT) image are given in 
lLopez et al.1 ([1995). All the images obtained at the 2.6 m 
Nordic Optical Telescope (NOT) were obtained using the 
same setup configuration (i.e. the Andalucia Faint Object 
Spectrograph and Camera, ALFOSC, and the [SII] filter 
centered on A = 6724 A and bandpass AA = 50 A) . More 
details on the acquisition and treatment of the NOT ob- 
servations can be found in lAnglada et al.l (|2007l ). Finally, 
the 2010 image was obtained at the 4.2 m William Her- 
schel Telescope (WHT), using the AC AM camera on the 
Cassegrain focus, giving a field of view of 8' with a spa- 
tial scale of ~ 0'.'25 pixel" 1 . A narrow-band [SII] filter, 
centered on A = 6727 A and bandpass AA = 48 A was 
used. In addition, a frame of a shorter exposure time of 
3600 s was acquired through another narrow-band filter, 
centered on A = 6645 A and bandpass AA = 50 A. This 
filter includes the nearby continuum, free of [SII] emission 
lines, useful to subtract the contribution of the contin- 
uum reflected emission close to the HH 30 jet source from 
the line emission image. For each epoch, several frames, 



Table 1 

[SII] images of the HH 30 jet/counterjet system 







Exp. time 




Epoch 


Telescope/Instrument 


(b) 


Reference 


1993 Dec 15 


INT 


14400 


1 


1998 Nov 08 


NOT/ALFOSC 


9000 


2 


1999 Nov 20 


NOT/ALFOSC 


14400 


2 


2006 Jan 19 


NOT/ALFOSC 


14400 


3 


2007 Nov 06 


NOT/ALFOSC 


10800 


3 


2010 Dec 01 


WHT/ AC AM 


12000 


3 



Ref erences. - - (1) lLopez et all (TMD (2) lAnelada etaT] 
B007t) (3) This work 



with a typical time exposure of 1800 s, were obtained 
to complete the total integration times, listed in Table 
[TJ of the final deep images. The individual frames were 
processed using the standard tasks of IRAlQ reduction 
package, which included bias subtraction and flat-field 
corrections, using sky flats. In order to correct for mis- 
alignments, all the individual frames of the same epoch 
were recentered using the position of field stars. Then, 
the frames were median-averaged to obtain a final deep 
image for each of the epochs listed in Table [TJ Images 
were not flux calibrated. 

The six final images were converted into a common ref- 
erence system and rebinned to the same pixel scale. The 
positions of seven field stars, common to all the frames, 
were used to register the images. The GEOMAP and 
GEOTRAN tasks of IRAF were applied to perform a 
linear transformation, with six free parameters that take 
into account translation, rotation and magnification be- 
tween different frames. All the transformed frames have 
a pixel size equal to that of the last epoch WHT image 
that was taken as the reference image. 

Astrometric calibration of the images transformed to 
the common reference system was made by using the 
coordinates of ten field stars, well distributed on the ob- 
served field. The coordinates were obtained after identi- 
fying the stars from the 2MASS All Sky Catalogue. The 
high optical extinction of the region made it impossible 
to find such a grid of reference stars from the USNO- 
B1.0 Catalogue. The typical rms of the transformation 
was < O'/l in both coordinates. The pixel size was found 
to be 0"2528 pixel" 1 . 

3. PROPER MOTIONS 

The images used for proper motion determination were 
rotat ed by an angle of 31 ? 6, the position angle of the jet 
axis (jAnglada et al.ll2"007fl , so that the y axis is along the 
jet axis and the x-axis is perpendicular, with values in- 
creasing from left to right. In order to improve the signal- 
to-noise ratio, the images were smoothed with a Gaussian 
with a FWHM of 3 pixels. The smoothed images were 
used for the farther away knots of the jet and counterjet 
(B to E, and J to N), while the full resolution images 
were used for the stronger, closer to the origin knots (A 
a nd Z). This nomencl ature for the knots is based on that 
of lLopez et all (|1995t) and lAnglada et al.l (|2007l) . Knots 
Z3 to Z6 of the counterjet, detected in our 2006, 2007, 

7 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 2 

Positions and Proper Motions of Knots in the HH 30 Jet 





x & 


2/ a 






fly 


e x b 








P.A. d 




Knot 


(") 


(") 


(" yr" 1 ) 




(" y^ 1 ) 


(") 


(") 


(km s~ 




(deg) 




Al 


-0.04 


2.69 


-0.005 ± 0.007 





.153 ± 0.010 


0.06 


0.08 


101.5 ± 


6.4 


1.9 ± 2 


.6 


A2 


0.08 


4.12 


0.007 ± 0.007 


0. 


150 ± 0.008 


0.08 


0.08 


99.5 ± 


5.1 


-2.6 ± 2 


.8 


A3 


0.09 


6.85 


0.003 ± 0.003 


0. 


179 ± 0.008 


0.03 


0.08 


118.7 ± 


5.1 


-1.0 ± 


.9 


Bl 


0.12 


11.04 


—0.005 ± 0.003 


0. 


097 ± 0.006 


0.03 


0.06 


64.5 ± 


3.9 


3.0 ± 1 


.9 


B2 


-0.10 


13.49 


-0.005 ± 0.002 


0. 


193 ± 0.007 


0.03 


0.10 


127.9 ± 


4.5 


1.5 ±0 


.7 


B3 


-0.25 


14.36 


0.008 ± 0.008 


0. 


190 ±0.010 


0.08 


0.10 


126.4 ± 


6.8 


-2.4 ±2 


.5 


C 


-0.30 


17.12 


-0.005 ± 0.003 


0. 


177 ±0.017 


0.04 


0.24 


117.3 ± 


11.4 


1.6 ±1 


.0 


Dl 


0.05 


21.19 


0.009 ± 0.007 


0. 


319 ± 0.028 


0.07 


0.29 


212.1 ± 


18.6 


-1.5 ±1 


.3 


D2 


0.06 


24.31 


-0.005 ±0.007 


0. 


382 ± 0.024 


0.10 


0.33 


253.8 ± 


15.8 


0.7 ±1 


.1 


D3 


0.09 


27.17 


-0.008 ± 0.005 


0. 


265 ± 0.098 


0.02 


0.34 


175.9 ± 


64.8 


1.8 ±1 


.2 


Dl 


-0.33 


29.64 




















El 


-0.60 


35.81 


-0.005 ± 0.002 


0. 


.149 ±0.007 


0.03 


0.10 


98.8 ± 


4.7 


2.0 ±0 


.8 


E2 


0.60 


41.05 


-0.001 ±0.009 


0. 


099 ±0.012 


0.12 


0.17 


65.9 ± 


7.9 


0.8 ±5 


.0 


E3b 


-0.55 


48.03 


-0.007 ± 0.009 


0. 


199 ± 0.004 


0.08 


0.04 


131.9 ± 


2.9 


1.9 ±2 


.5 


E4 


-0.46 


52.04 


-0.020 ±0.007 


0. 


.143 ±0.003 


0.10 


0.05 


95.8 ± 


2.3 


7.8 ±2 


.7 



a Positions in the 2010 image, except for knot B3, for which the position is that of the 2007 
image. The y axis is along the jet, at a position angle (eastwards from north) of 31?6. The 
(0, 0) position is that of the brightest knot, AO. 

b Rms residual of the knot positions in the proper motion fit, e = ay/1 — r 2 , where a is the 

standard deviation and r the correlation coefficient. 

c Proper motion velocity, assuming a distance of 140 pc. 

d Position angle with respect to the y direction. 



Table 3 

Positions and Proper Motions of Knots in the HH 30 Counterjet 





X 3, 






My 




% b 




P.A. d 




Knot 


(") 


(") 


(" yr- 1 ) 


(" yr- 1 ) 


(") 


(") 


(km s _1 ) 


(deg) 




Zl 


0.05 


-5.61 


0.000 ±0.010 


-0.278 ±0.012 


0.08 


0.13 


184.8 ± 8.3 


-179.9 ± 


2.1 


Z2 


-0.05 


-9.55 


0.011 ±0.001 


-0.462 ± 0.003 


0.01 


0.03 


307.0 ± 2.2 


-178.7 ± 


0.1 


Z3 


0.16 


-11.71 


0.065 ±0.020 


-0.223 ±0.026 


0.07 


0.09 


154.4 ± 17.1 


-163.7 ± 


5.1 


Z4 


-0.23 


-13.43 


-0.022 ±0.022 


-0.166 ±0.055 


0.08 


0.19 


111.0 ±36.0 


172.5 ± 


8.0 


Z5a 


0.06 


-17.41 


0.046 ± 0.046 


-0.227 ±0.003 


0.16 


0.01 


153.8 ± 6.4 


-168.6 ± 


11.3 


Z5b 


-0.10 


-18.26 


-0.027 ±0.006 


-0.316 ±0.003 


0.02 


0.01 


210.7 ± 1.9 


175.1 ± 


1.0 


Z6 


0.16 


-23.65 


0.012 ±0.016 


-0.608 ± 0.084 


0.05 


0.29 


403.3 ±55.5 


-178.9 ± 


1.5 


J 


0.32 


-78.30 


0.029 ±0.013 


0.066 ± 0.060 


0.17 


0.78 


48.0 ±36.9 


-23.6 ± 


21.5 


Kl 


1.56 


-94.35 


0.064 ±0.016 


-0.174 ±0.057 


0.20 


0.73 


122.9 ±35.5 


-159.8 ± 


7.6 


K2 


0.89 


-99.02 


0.036 ±0.012 


-0.148 ±0.032 


0.16 


0.42 


101.3 ±21.0 


-166.4 ± 


5.4 


L 


4.09 


-149.00 


0.004 ±0.012 


-0.249 ±0.025 


0.16 


0.32 


165.0 ±16.6 


-179.2 ± 


2.8 


M 


7.23 


-176.61 


0.011 ±0.008 


-0.325 ±0.015 


0.10 


0.20 


215.9 ±10.2 


-178.0 ± 


1.4 


N 


11.59 


-191.00 


0.076 ± 0.003 


-0.323 ±0.015 


0.04 


0.19 


220.1 ± 9.5 


-166.7 ± 


0.8 



a Positions in the 2010 image. The y axis is along the jet, at a position angle of 31?6. The (0,0) 
position is that of the brightest knot, AO. 

b Rms residual of the knot positions in the proper motion fit, e = ay/1 — r 2 , where a is the standard 
deviation and r the correlation coefficient. 
c Proper motion velocity, assuming a distance of 140 pc. 
d Position angle with respect to the y direction. 

and 2010 observations, have been identified for the first 
time. Positions of the knots for each epoch were mea- 
sured with respect to the position of the brightest knot 
AO, whose position was set to (x,y) — (0,0). The po- 
sition of the knots was determined from a parabolic fit 
to the intensity of the 5x5 pixels centered on the pixel 
with peak emission. Then, the proper motion in x and y 
directions of each knot, \x x and fi y , was determined from 
a linear regression fit to their positions in the different 
epoch images. 

The 1993 image has poorer seeing than the others, and 
could not be used for knots close to the HH 30 star, i.e. 
knots A and Z. The 1989 and 1999 images do not include 
the counterjet region for distances greater than 5'.'6, and 
could not be used neither for the counterjet knots Z3 to 



Z6, nor for J to N. The results are given in Tables [5] and [31 
and shown in Fig. [JJ where the knot positions and proper 
motions are shown superimposed on the smoothed 2010 
image. 

The errors assigned to the proper motions, in both x 
and y, are the formal errors of the slope of the linear 
regression fits for each knot. These errors are shown as 
ellipses at the end of the arrows shown in Fig. [TJ The 
errors appear to be rather small as a consequence of the 
wide span of time used for calculating the proper mo- 
tions, and that 6 epochs were used for most knots. The 
quality of the linear regression fits for each knot is indi- 
cated by the small values of the residuals in x, e x , and 
y, e y , (see Tables [3] and [3]) and is illustrated in Fig. [5J 
where the measured and fitted values of the y position 



X offset (arcsec) 



X offset (arcsec) 



X offset (arcsec) 



Figure 1. 2010 [SII] narrow-band image of the jet and counterjet of HH 30, smoothed with a 3-pixel Gaussian, showing the proper motion 
velocities of the knots A to E (left); knots Z (center); and knots J to N (right). The y axis is along the jet axis, at a position angle of 31?6. 



as a function of time for each knot, are shown. In all the 
cases, the y positions are well fitted by the proper motion 
velocity obtained, showing that the motion of the knots 
is ballistic. 

All the knots show proper motions roughly in the direc- 
tion of the jet axis, with the exception of knot J of the 
counterjet, which has been measured to be nearly sta- 
tionary. This knot is located inside a reflection nebulosity 
southwest of the HH 30 star (spreading from y ~ —25" 
to —100"; see Fig. [TJ, and is most probably not a coun- 
terjet knot, but a feature of the reflection nebulosity. We 
will not consider this knot in the following. 

It is worth noting that the proper motions ji y ob- 
tained here are signi ficantly lower than those obtained in 
lAnglada et al.l ()2007t ) from two observations with nearly 
one year of interval, in 1998 and 1999. On the average, 
the proper motion velocities for kn ots A to E of the jet are 
O'.'ll yr" 1 (or 75 km s -1 ) higher in lAnglada et al.l (|2007h 
than the present values. For the counterjet, the absolute 
value of the proper motion velocities for knots Zl and 
Z2 are C/13 yr _1 (or 85 km s _1 ) lower in lAnglada et alJ 
(2007) than the present values. Since the proper mo- 
tions are measured with respect to the brightest knot 
AO, the discrepancy can be explained by a shift in the 



position of AO along the axis of the jet of Ay ~ — 0"12 
in the 1999 image with respect to that of 1998. Such a 
shift in position can be attributed to the observed vari- 
ations i n the brightness of the diffuse light of the HH 
30 st ar (|Watson fc Stapelfeldtl l2007t iDuran-Roias et alJ 
2009). However, this systematic error in the determina- 
tion of proper motions is not expected to affect measure- 
ments that span over a large number of years, like those 
presented in the present work. 

In Fig. [3] we plot the proper motion velocity in the y 
direction, in absolute value, | function of distance 

y. As can be seen, the jet and counterjet velocities are 
quite different. For the jet knots, most have velocities 
~ 100 km s _1 , except knots D, with distances from the 
HH 30 star between 20" and 30", which have a higher 
velocity, ~ 240 km s" 1 . For the counterjet, the velocities 
are more irregular. There are knots with \v y \ ~ 100 km 
s- 1 (Z4), with \v y \ ~ 180 km s" 1 (Zl, Z3, Z5a, Z5b), 
and with \v y \ > 300 km s _1 (Z2, Z6). On the average, 
the proper motion velocity of the jet knots is 115 km 
s _1 , while for the counterjet knots the average is 193 km 
s _1 . The difference is also found in the velocity disper- 
sion, which is 54 km s _1 for the jet, and 82 km s _1 for 
the counterjet. Thus, the counterjet knots have, on aver- 
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Table 4 

Launch properties of knots in the 
HH 30 jet and counterjet 





f a 
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Figure 2. Diagram showing the y position of the knots of the jet and counterjet, as a function of time: knots A to E (ie/i); knots Z 
(center); and knots J to N {right). The red lines show the linear fits used to derive the proper motion velocity in the y direction. The 
positions of 1993 have not been used for the proper motion determination of knots Z. The 1998 and 1999 images only reach up to y = — 5'.' 6 
from the HH 30 star in the counterjet. 

age, higher velocities than the jet knots, and show more 
variation from knot to knot. 

Let us now consider the pair of knots of the counter- 
jet Zl and Z2. The velocity of Z2 (307 ± 2 km s" 1 ) is 
much higher than that of Zl (185 ± 8 km s" 1 ). This 
pair of knots, in spite of being very close to the HH 30 
star, have well determined proper motions, since they 
have been observed at five epochs, spanning 17 years. 
In addition, the proper motions obtained by fitting only 
the new observations in 2006, 2007, and 2010, coincide 
with those obtained when the old observations of 1998 
and 1999 are included in the fit. The positions of the 
two knots Zl and Z2, when extrapolated back in time, 
show that they were launched simultaneously, at the end 
of 1989 (see Fig. [2]). A similar case is found for the pair 
of knots Z5 and Z6 of the counterjet, but in this case the 
proper motion velocities are not so well constrained as in 
the case of Zl and Z2. Thus, it appears that the coun- 
terjet of HH 30 is able to launch simultaneously pairs of 
knots with very different velocities. 

In order to investigate the launching properties of all 
the knots, we calculated for each knot its kinematic age, 
tkin ~ y/Hy, and the position at launch, xi c h — x— /^ikin- 
The results obtained are shown in Table [4] In Fig. [4] we 
show the x position at launch, as a function of dis- 
tance to the origin, \y\. As can be seen, all knots, except 
the more distant knots, have values of xi c h consistent 
with zero. This is an indication that these knots are bal- 
listic, moving along a straight line after their launching. 
For the more distant knots, Kl, K2, and N, x\ c h — —30", 
indicating that they deviate from a straight line after 
their launching. As we will see in Q the more distant 
knots trace the large-scale "C" -shape bending of the jet 
and counterjet. 
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Figure 3. Proper motion velocity in the y direction (absolute 
value), \v y \, versus distance, y, for the knots in the HH 30 jet 
(blue) and counterjet (red). The stepped lines indicate the veloc- 
ities adopted for the knots of the jet (blue) and counterjet (red). 
These velocities are the error-weighted averages of the velocity of 
the knots in each step: Vj = 98 km s -1 (A, B, C, E, Z4); Vj = 181 



km s" 1 (Zl, Z3, Z5); Vj 
(7,2, Z6). 



■ 234 km s" 1 (D); and Vj = 307 km s" 



4. JET BENDING 

4.1. Bending by a plane-parallel side wind 

The bending of the HH 30 je t/counterjet, with a clear 
"C" shape, was rep orted by lLopez et al.l (|1995l ) and 
lAnglada et al.l (|2007l) . Now, with the observations of the 
counterjet presented in this work, we know the shape of 
the jet/counterjet for a length of ~ 500", corresponding 
to ~ 0.35 pc. 

The bending of the HH 30 jet is suggestive of being 
entrained toward the north-west (increasing x). This 
could be due to a proper motion of the HH 30 star to- 
ward the south - east w ith respect to the ambient medium. 
ICanto fcT laga (1995) estimate that a relative velocity 
between the HH 30 star and the surrounding environ- 
ment of ~ 2 km s _1 is sufficient for explaining the ob- 
served bending of the HH 30 jet/counterjet. Such a rel- 
ative velocity could be the result of a proper motion of 
the HH 30 star of ~ 0'.'003 yr" 1 , clearly below our ac- 
curacy in the proper motion measurements, of the order 
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Figure 4. Position in the x axis at the epoch of launch assuming 
ballistic motions, Xi ch = x — y fi x //i y , versus \y\ for the knots in 
the HH 30 jet (blue) and counterjet (red). Top: All knots. The 
dashed line shows the region that appears enlarged in the bottom 
panel. Bottom: Knots close to the origin. 

of ~ O'.'l yr^ 1 . An alternative is that the deflection is 
caused by the powerful bipolar outflow whose exciting 
source, L1551-IRS5, is at a distance of ~ 4' southwards 
of HH 30. The red lobe of the L1551-IRS5 molecular 
outflow is found to the south-east of the HH 30 star, and 
could be responsible of entraining the surrounding envi- 
ronment with a velocity of ~ 2 km s^ 1 and deflecting the 
jet and counterjet of HH 30. 

A different explanation for the bending of the HH 30 
jet and counterjet is a deflection by an isotropic stellar 
wind blowing the HH 30 jet from the south-east. In the 
following we will examine this possibility. 

4.2. Bending by an isotropic stellar wind 

We will use the model o f jet deflection by a n isotropic 
stellar wind, developed by iRaga et al.l j2009a) , to fit the 
shape of the HH 30 jet/counterjet. The path of the 
jet/counterjet in a reference system centered on the wind 
source (see Fig. [5]), is given in polar coordinates by 



<<t>) 



cos A (0/A) ' 



A 2 < 



(i) 



where r s is the distance from the wind source to the 
"stagnation point" (i.e. the point closest to the wind 
source) of the jet/counterjet system, and 



A 



< A < 1, 



with e given by 



e = 2 



1/2 



(2) 



(3) 
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Figure 5. Reference systems used for studying the jet and coun- 
teract bending. The jet path is indicated by the thick continuum 
line. The jet source, J, is indicated with a small circle, while the 
wind source, W, is indicated with a small square. The "stagnation 
point" is S. The distance from the side-wind source to S is r B , and 
to J is v w . The angle ct w is the angle between the x f and y axes, 
and <f>s is the angle with vertex at W, subtended by the jet source 
J and the stagnation point S. The reference system (x',y') is ro- 
tated an angle <j> s with respect to (x",y"). The reference system 
(x, y) is shifted a distance r m and rotated an angle ct w — ir/2 with 
respect to (x',y'). The two dashed lines indicate the asymptotes 
to the jet and counterjet. The total deflection of the jet, ipj , is the 
angle between the two asymptotes. 

where Mj, Vj and Cj are the mass-loss rate, velocity and 
isothermal sound speed of the jet, and M w , v w are the 
mass-loss rate and velocity of the wind. The parameter 
A can be expressed in terms of the total deflection angle 
of the jet/counterjet path, tpj (see Fig. [5]), 



X = l-tl. 

IT 



(4) 



Let us call r w the distance from the jet source to the 
wind source, and a w the angle between the jet axis {y 
axis) and the direction from t he wind source to th e jet 
source (a;' axis). As shown by iRaga et al.l (|2009aD . the 
angle <j) s , with vertex at the wind source and subtended 
by the jet source and the stagnation point (see Fig. [5]), 
is given by. 

77 

(5) 



K a w - 1) 



The polar equation of the jet shape in the reference sys- 
tem (x\ y') is given by 



= 



sin a v 



sin(a w — 4>' /X) 



X(a w — 7r) < 4>' < Xa v 



(6) 

Therefore, the coordinates of the jet shape in the refer- 
ence system (x,y), centered on the jet source and with 
the y axis along the tangent to jet/counterjet path, are 
given by 



r(cj>') sin <fi' 



cos a v 



x = [r (</>') cos 4>' — r w \ sin a u 
y— [r ((/>') cos0' — r w ] cosa^, + r(<p') sin<p' sina^. (7) 

The jet shape described by Eqs. [6] and [7] depends on 
three parameters: the distance to the wind source, r w ; 




X offset (arcsec) 

Figure 6. Bending of the HH 30 jet and counterjet. Right: The 
small circles mark the positions of the knots of the HH 30 jet 
and counterjet, with the farther knots labeled. The continuum 
line sho ws the best fit to the knot positions using the R aga et alj 
( 2009a) model of jet bending by an isotropic stellar side wind. The 
x scale is amplified a factor of ~ 10. Left: Best fit superposed on a 
[SII] i mage of the jet and c ounterjet. The upper half-image is from 
1999 (Anglada et al. 2007) and the lower half-image is from 2010 
(present work). The x and y scales are the same. 

the angle between the jet axis and the direction to the 
wind source, a w ; and the total deflection angle, ipj (or 
A, related through Eq. EJ). 

A least-squares fit of the three free parameters of the 
model w as performed to the po sitions of 43 knots, NA to 
NH from lAnglada et al.l (|2007D . knots Al to 12 of the jet, 
and Zl to N of the counterjet (present work), spanning 
a total length along the jet axis of ~ 500". 

The best fit was obtained for r w = 173"±15" (24200± 
2100 AU in projection), a w = 66° ± 4°, ip = 14?8 ± 0?7 
(corresponding to A = 0.918 ± 0.004), with an rms fit 
residual in x of eg t = 2'.'22 (see Fig. [HJ) - The errors quoted 
were estimated as the increment in each parameter that 
increases the rms fit residual a factor of (1 + "X 2 3/n) 1 / 2 , 
where n is the number of points fitted, and x 2 3 is the 
value of x 2 f° r 3 degrees of freedom (the number of 
free parameters) and 68% significance (l-er uncertainty), 
X 2 3 = 3.53 (jLampton et al.fl l976). For these values of the 
parameters we obtain that the stagnation point is located 
at (x, y) = (+0'.'8, -64") (5 position in Fig.©. The wind 
source is located at a position (x, y) = (—158", —70") (W 
position in Fig. [7]). Within the 1-a uncertainty ellipse of 
the W position we found a reddened star, identified in the 
2MASS catalog (J04314418+181047), with coordinates 
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Figure 7. Bending of the HH 30 jet and counterjet by an 
isotrop ic stellar wind. The dotted line shows the best fit using 
the Raea ct al. (2009a) model. The positions of the HH 30 star 
(J), wind source (W) and stagnation point (5) are indicated. The 
wind source is located at an offset position (—158", —70"). The 
ellipse shows the 1-cr uncertainty (68% confidence) in the position 
of W. The stagnation point is located at (+0"8, -64"). 

a(J2000) = 04 h 31 m 44!180, <5(J2000) = +18°10'14"78, 
whose colors (J — H = 1.33, H — K = 0.81) corre- 
spond to those of a Classical T Tauri star (CTTS) in 
the (J — H),(H — K) diagrams (|Lada fc Adam s 1992; 
IRoman-Lopesll2009l) . 

The total deflection of the jet, ip = 14?8±0?7, implies 
a value of the parameter e — 11.2 ± 0.6. As stated in Eq. 
02 e is related to the physical properties of the jet and the 
wind source. Taking a jet velocity Vj — 100 km s , and 
an isothermal sound speed of the jet Cj — 10 km s -1 , we 
obtain 



M jVj 



2 v« 



= 0.31 ±0.02. 



(8) 



The mass-loss rate of the HH 30 jet and counterjet 
has been estimated t o be Mj ~ 2.6 x 10 -9 Mq yr _1 
(jBacciotti et all I1999D . For a jet velocity Vj = 100 
km s _1 , the HH 30 momentum rate is MjVj ~ 2.6 x 10~ 7 
Mq yr _1 km s" 1 . Thus, the wind source needed to de- 
flect the jet has to have a momentum rate 

M w v w ~ (8.3 ± 0.5) x 10~ 7 Mq yr^ 1 km s _1 . (9) 

Mass-loss rates for CTTS in Taurus are of the order of 
M w ~ 10 _10 -10~ 8 Mq yr^ 1 , and up to one order of mag- 



nitude higher for continuum CTTS (|Hartmannl 120051 ). 
Typical CTTS wind velocities are of the order of 100 
km s _1 . Thus, for a CTTS we can expect 



M w v w ~ 10~ 8 -10 -6 M Q yr" 1 km s 



-l 



(10) 



This range of momentum rates encompasses the value 
needed to cause the observed deflection of the jet, given 
by Eq. 9. Thus, the star J04314418+181047 could be 
the wind source responsible for the deflection of the HH 
30 jet, provided its mass-loss rate falls in the higher end 
of values found for CTTS. 

If a CTTS can produce the observed bending of HH 30, 
this would imply that a vast majority of jets should be 
bent, since low-mass protostars are abundant. However, 
the bending of a jet is not easy to detect since it is only 
noticeable when the jet is imaged over a long length. In 
the case of HH 30, the images span a length of 500", or 
0.35 pc. 

5. JET WIGGLING 

Following the work of lAnglada et al.l (|2007l ), we con- 
sider that the wiggling of the HH 30 jet and counterjet 
is a consequence of the presence of a companion star 
to the jet source. lAnglada et al.l (l2007h. u sing the for- 
mulation given byUVlasciadri & Raga (2002), analyze the 
extreme cases where the dominant effect is either the 
orbital motion of the jet source in a binary system, or 
the precession of the ejection axis of the jet because of 
tidal interactions between the disk where the jet origi- 
na tes and a companio n star. Here we develop the work 
of iRaga et al.l ()2009bl ). and we consider a physical sys- 
tem, in which both the orbital motion and precession 
can be present simultaneously, to fit the HH 30 jet and 
counterjet wiggling shape. 

5.1. Orbital and precession periods 

We consider a binary system with a circular orbit, be- 
ing uij the mass of the jet source, m c the mass of the 
companion, and m = rrij + m c the total mass of the sys- 
tem. We will call /i c the mass of the companion relative 
to the total mass, so that 



rrij = (1 - jit c ) m, 
m c = fx c m. 



(11) 

Let a be the binary separation (i.e. the radius of the 
relative orbit). Therefore, the orbital radius of the jet 
source with respect to the binary's center of mass (i.e. 
the radius of the jet source absolute orbit) is 

r = [i c a, (12) 

and the orbital velocity of the jet source is given by 

v = , (13) 

To 

where r Q is the orbital period. The total mass of the 
binary system is related to r , and /i c , r , or a, through 
Kepler's third law, 
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yr 



(14) 

Let us consider that the disk of the jet source is tilted 
an angle j5 with respect to the orbital plane, and that it is 
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Figure 8. Reference systems used for studying the jet and coun- 
terjet wiggling. The jet axis has an inclination i with respect to 
the plane of the sky. The orbit is contained in the (x' ,y') plane. 
The orbital phase angle, O , is measured counter-clockwise from 
the x' axis. 

precessing with a period t p . An approximate expression 
relating the or bital and preces sion periods can be derived 
from Eq. 24 of lTerqueml (|1998D , valid for a disk precessing 
as a rigid body, by assuming that the disk surface density 
is uniform and that the rotation is Keplerian, 



15 



32 (1-Mc) 1/2 



(T 3 / 2 COS/ 



(15) 



where a = r^/a is the ratio of disk radius to binary 
separation. Since it is expected that the size of the disk 
is truncated by tidal interaction with the companion star 
in such a way that 1/4 < a < 1/2 (|Terquemlll998D . we 
will adopt a value of a = 1/3. With this value of a, Eq. 
15 gives: t /t p = 0.09(^/^1 — Mc ) cos ft. Note that, 
except for values of fi c near to 1, the precession period is 
much longer than the orbital period. 

5.2. Velocity and position of the knots as a function of 

time 

For describing the jet we use a (linear) coordinate sys- 
tem (x',y',z'), with origin in the jet source, and where 
(x 1 , z) is in the orbital plane, the cc'-axis is the intersec- 
tion of the orbital plane with the plane of the sky, and 
the y'-axis coincides with the orbital axis, at an inclina- 
tion angle i with respect to the plane of the sky (pointing 
away from the observer). The coordinate system is illus- 
trated in Fig. ISJ 

As in Sj3] for describing the observations we use an 
(angular) coordinate system (x, y, z), where (x, y) is in 
the plane of the sky (the y-axis in the jet direction and 
the x-axis perpendicular with values increasing from left 
to right), and the z-axis points away from the observer 
(see Fig. [5]). The transformations between both systems 
of coordinates are 



x' 
D' 



V 



y cos z — z sin i 
D ' 



(16) 



where D is the distance from the source to the observer, 
which accounts for the linear to angular transformation. 

Let us assume that the velocity of the jet has a com- 
ponent perpendicular to the disk plane, Vj, and a com- 
ponent in the orbital plane, caused by the orbital motion 



of the jet source, v Q . The components of the velocity of 
a jet parcel, or knot, are 



v x ' — —Vo sin 9 D ± Vj sin ft cos 9 p , 
v y i = ±Vj cos ft, 

V, 



(17) 



= v Q cos 9 ± Vj sin ft sin 6 p , 



with the sign of the Vj term being positive for y' > 0, 
and negative for y 1 < 0. The angles 8 Q and 8 p are the 
orbital and precession phase angles, measured from the 
x' axis (perpendicular to the line of sight), at the epoch 
of launch, t\ c h, when the jet parcel was ejected, 



2tt 



(*lch — to) + <Po, 



2tt 



(tich — to) + (p p , 



(18) 



where tp a , ip p (between and 2n) are the orbital and pre- 
cession phase angles at an arbitrary epoch of reference, 
to- Note that we assume that the precession is retrograde 
(see t he discussion on ret rograde precession in binary sys- 
tems; Montgomery 2009), i.e. the orbital phase angle 9 
increases with time, while the precession phase angle 9 p 
decreases with time. The expression of Eq. [18] corre- 
sponds to orbital phase angles increasing from x' to z' 
(see Fig. [5]), i.e. a counter-clockwise rotation as seen from 
positive values of y' (and clockwise for precession). In 
case of clockwise orbital rotation (and counter-clockwise 
precession), the orbital phase angle would decrease from 
x' to z' , and the velocity v z > (and the z' coordinate, see 
below), would have the opposite sign. The difference is 
only noticeable when the inclination angle i is large, since 
z' appears multiplied by sin i in the transformation given 
by Eq. [TH 

The velocity of a jet parcel, or a knot, in the plane of 
the sky can be found from Eq. [17] by using a coordinate 
transformation similar to that of Eq. [T6l 

v x = — v Q sin o ± Vj sin /3 cos B p , 

v y = —v cos 9 a sin i ± Vj(cos ft cos i — sin ft sin# p sin ($9) 

The trajectory of this jet parcel, or knot, as a function 
of time is a straight line given by 

x = r a cos o + (~v a sin O ± Vj sin ft cos 9 p ) (t — ti c h), 
y' = ±v ]C osft(t-t lch ), (20) 
z = r a sin o + (v a cos 9 Q ± Vj sin ft sin 6 p ) (t — ii c h). 

Note that the position of a knot at the epoch of its ejec- 
tion (t = ti c h) is in the orbital plane (y' = 0). By using 
the coordinate transformation of Eq. [16] we obtain that 
the position of the knot in the plane of the sky, as a 
function of time, is given by 

x = — [r Q cos 8 + (—v sin 9 a ± Vj sin ft cos 9 P ) (t — ti c h)] , 



y = — [— r a sin 9 a sin i + {—v Q cos 9 Q sin i 

± Vj [cos ft cos i — sin ft sin 9 P sin i]) (t — tich)] ■ 



(21) 



5.3. Constant velocity jet: shape of the jet 

If the jet velocity Vj is constant, all the knots have the 
same velocity, and the jet at a given epoch, t = t b s , has 
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a shape that can be obtained by eliminating the time de- 
pendence of Eq. [20l substituting t — t\ c h = t b s — Uch = 
\y'\/(vj cos/?) in the equations. The orbital and preces- 
sion phase angles can be expressed as 



9 - 2n U 

' o \ ^obs 



2tt 



to 



V 



Vj COS /3 

\y'\ 



■<Po, 



Vj cos f3 



(22) 



For simplicity, we can introduce the parameter a (which 
plays a role similar to the angle of precession /3), 



tana 



Vj cos j3 ' 



(23) 



so that v (t b s — ti c h) = \y'\ tana. The angle a is the half- 
opening angle of the jet in the (x',y') plane, containing 
the jet axis, in the absence of precession. With these 
changes, the shape of the jet is given by 

x = r Q cos 9 Q — \y\ tan a sin 9 Q + y' tan /3 cos 9 P , 

z' = r Q sin 9 a + \y'\ tan a cos 9 + y' tan (3 sin9 p . (24) 

The shape of the jet in the plane of the sky (x, y) is 
obtained from the last equation, using the transformation 
of Eq. US 

5.4. Variable velocity jet 

The proper motions measured for the HH 30 knots (<j3]) 
show that the ejection velocity is significantly different 
for some knots. Let us consider that Vj (i.e. the ejection 
velocity component perpendicular to the disk plane) , in- 
stead of being the same for the full jet, can be different 
for different parts of the jet. In this case, the jet has 
a shape composed of different parts moving at different 
velocities. Let us estimate the x position of each knot 
predicted by a model with given values of the parame- 
ters of the orbital motion (v , t , fo), and of those of the 
precession (/3, t p , tp p ). 

For each knot we know, from our observations, its po- 
sition, y, and proper motion, fj, y , along the jet axis. We 
can approximate the knot kinematic age by t^ ln ~ y/l^y 
This approximation does not take into account the shift 
in position of the jet source caused by its orbital mo- 
tion, so it holds for y 3> r a /D, or for small inclination 
angles i. The epoch of launch of the knot is given by 
t\ch = tohs — tkinj and we can estimate the orbital and 
precession phase angles at the epoch of launch, given by 
Eq. [THJ Once the orbital and precession phase angles are 
known, the ejection velocity component perpendicular to 
the disk plane for this knot, Vj, can be calculated from 
Eq. OH 



Dfly 



v cos t) sin % 



cos /3 cos i — sin /3 sin 9 P sin i 



(25) 



Therefore, the position of the knot in the plane of the sky 
at the epoch of observation (t — t b s ) can be obtained 
from Eq. [2U for t - t ich = t kin , 

x~ — [r cos9 Q + {—v sin0 o + Vj sin (3 cos 9 p ) t kin ], 



= y - 



Vj [cos (3 cos i — sin /3 sin 9 p sin i] ) i k ; 
D 



■ sm v sm i . 



The first equation gives the x coordinate of the knot, and 
the second equation shows that the y coordinate of the 
knot is derived consistently by the model, for y 3> r a /D, 
or for small inclination angles i. 

5.5. Model parameters and limiting cases 

In general, if the inclination angle i is known, and both 
orbital motion and precession can be fitted to the data, 
the model depends on six parameters, which can be cho- 
sen to be the mass function, mp? c (see below), the orbital 
period, r Q , the orbital phase angle at epoch t$, ip a , the 
precession angle, /3, the precession period, r p , and the 
precession phase angle at epoch to, if p - The rest of pa- 
rameters of the system can be obtained from these: /i c 
is derived from Eq. 1151 rrij = (1 — /i c )m, and m c = fi c m 
(Eq. OH); r Q and a are derived from Eq. 1141 and v a is 
derived from r D and r Q (Eq. [T3|) . 

The "pure" orbital model of § 4.1.1 of I Anglada et all 
(|2007| ) corresponds to j3 — in the previous equations 
(t p , and tp p are undefined). In the case of a constant 
velocity jet, the shape of the jet (Eq. |2~4"|) has a plane 
symmetry ("C" symmetry) with respect to the x axis. 
For such a "pure" orbital model, a (Eq. |23|) is the half- 
opening angle of the jet cone, and the orbital period r D 
is determined by the jet velocity and the wavelength of 
the jet wiggles in the plane of the sky. However the total 
mass of the system, m, and the mass of the companion 
relative to the total mass, /i c , cannot be determined in- 
dependently, since the model depends only on the mass 
function m/j^. (Eq. [T4|) . Thus, the "pure" orbital model 
depends on three parameters, which can be chosen to be 
the mass function, irifx^., the orbital period, t q , and the 
orbital phase angle at epoch to, tp a . This will be also 
the case of the general model when precession is not rel- 
evant in shaping the observed jet, for instance when the 
precession angle /3 is small compared to a, or when the 
wavelength of the wiggles caused by precession is longer 
that the length of the observed jet. 

The "pure" prec ession model of § 4.1.2 of 
lAnglada et al.1 (|2007ft is obtained from the previous 
equations in the limiting case /i c — > 0. For this limit 
the orbital radius r Q , the orbital velocity v a , and a 
tend to zero, while the binary separation a and orbital 
period t d remain finite. In this case the orbital terms 
in Eqs. |2"41 and l2l)l vanish, and, in the case of a constant 
velocity jet, the shape of the jet has a point symmetry 
("S" symmetry) with respect to the origin. However, 
for a real system (/i c > 0), there is always orbital 
motion, even if precession is dominant. For such a 
precession-dominant model (a -C ft), the precession 
angle (3 is the half-opening angle of the jet cone, and the 
precession period t p is determined by the jet velocity 
and the wavelength of the jet wiggles in the plane of 
the sky. Thus, three parameters, the precession angle, 
/3, the precession period, t p , and the precession phase 
angle at epoch to, (p p , can be obtained from the fit. 
Then, for each value of /i c in the interval < fj, c < 1, 
an orbital period r Q is obtained from Eq. [15] and, since 
a <C P, the following upper limits are found for the 



The Counterjet of HH 30 



11 



binary separation, a, (Eqs. IT21 [TBI [151 |2"3"|) 
0-09 1 . „ 

a <<: (i- Me )i/2 t p ^ sm p cos & ( 27 ) 

and for the mass of the jet source, m,j, (Eqs. [TTJ [TJ] [15j 

EH) 

m j < 7TT ~ " ^ 7 T P v i sin3 Z 3 cos P- ( 28 ) 

5.6. Fitting procedure 

For the fitting process we consider as known the in- 
clination angle of the jet axis with respect to the plane 
of the sky, i. As di scussed in lLopez et al.l (|1996l) and 
Angla da et al.l (|2007T > the inclination angle is small, and 
so, the model depends weakly on i. We have taken i = 5°. 
Another parameter considered as known is the offset in 
the y direction of the HH 30 star position with respect to 
the intensity peak of HH 30 (knot AO). W e have taken 
— 0'.'51, the same value as in lAnglada et all (|2007t) . 

The model of variable velocity jet was used to fit the 
x positions of the knots A to E of the jet and the knots 
Zl to Z6 of the counterjet in the 2010 image (Tables [5] 
andEJ. The knot D4 of the jet was excluded from the fit 
because it was identified only in the 2010 image and its 
proper motion could not be measured. 

The fitting strategy was to sample the six-dimensional 
parameter space, defined by the parameters m/ijj, r D , 
/?, (fi and ip p . Several sampling methods were tested, 
i.e. r egular grid, random, and Halton sequence (jHaltonl 
Il964f ). We discarded the regular grid method because 
some parameters (for instance, r Q ) can have evenly 
spaced values that fit the data. We adopted a Halton 
quasi-random sequence because it samples the space pa- 
rameter more evenly than a purely random sequence, and 
the convergence of the fitting procedure to the minimum 
of the rms residual is faster. 

Several runs of the fitting procedure were performed 
using different number of sample points (typically be- 
tween 10 6 and 10 7 ), and different ranges of the parame- 
ters. The "pure" orbital model, depending three param- 
eters, mp? c , t , and <£> was also fitted, by adopting a fixed 
value (3 — 0. 

Once a minimum of the rms fit residual was found, the 
uncertainty in the parameters fitted was found as the 
increment of each of the six (or three) parameters of the 
fit necessary to increase the rms fit residual a factor of 
(1 + x 2 p / n ) 1 ^ 2 : where n is the number of knots whose 
position was fitted, p is the number of parameters fitted, 
and x 2 p is the value of \ 2 f° r P degrees of freedom (the 
number of free parameters) and 68% significance (1-cr 
uncer tainty) , i.e. x 2 3 = 3.53, x 2 q = 7-04 (jLampton et al.l 
[19761) . 

5.7. Results of the model fitting 

Let us first consider a system where precession and or- 
bital motion are present in the jet, and both contribute 
to determine the shape of the jet. This condition means 
that the half-opening angles due to orbital motion, a, 
and to precession, /3, are of the same order. For the 
range of masses considered for the HH 30 jet source and 
the companion (0.06-1 M Q ; see below), the orbital pe- 
riod is substantially shorter than the precession period. 




m/JT 1 



Figure 9. Rms fit residual of the "pure" orbital model fitted to 
the HH 30 knots of the jet and counterjet, as a function of mass 
function mfj!^ and orbital period t . The cross marks the position of 
the best fit, for which the rms is 0'.'104. The blue contours indicate 
the confidence levels 68% (1-ct uncertainty), 90%, and 99%, while 
the red contour indicates the rms x position of the knots with 
respect to the jet axis, 0'.'272. 

Thus, the jet shapes obtained from the model present a 
short-scale wiggling corresponding to the orbital period 
T , and a large-scale wiggling corresponding to the pre- 
cession period t p (see, for instance, Fig. 3 of iRaga et al.l 
2009b). Although such a kind of model can reproduce 
the observed x positions of the limited number of knots 
included in the fitting, it predicts a small-scale wiggling 
shape between these knots, which is not observed in the 
HH 30 jet images. Thus, the model that fits the data 
has to be either precession-dominant (a <C or orbital- 
motion-dominant < a). 

The precession-dominant model can be excluded from 
simple considerations. Since the observed wiggles of 
the jet in the plane of the sky have a wavelength of 
~ 16" (jAnglada et all 120071 see also Fig. 14), a jet ve- 
locity of Vj ~ 100 km s _1 would result in t p ~ 100 
yr, assuming a precession-dominant scenario. Since the 
observed half-openin g angle of the jet cone is ~ 1?5 
(Angl ada et al.l [20071) . and a <§; (3, the precession an- 
gle should be (3 ~ 1?5. We know that the mass of the 
jet source cannot be higher than ~ 1 M®. Higher values 
of the mass are unlikely since th e bolometric lumi nosity 
of the system is less than 1 L Q (|Cotera et al.|[200ll) . It 
can neither be lower than the mass of a brown dwarf, 
~ 0.06 Mq. Using these constraints and Eqs.[27l[28l we 
can derive an upper limit for the binary separation for 
the allowed range of masses of the jet source. The value 
obtained for the binary separation in all cases is a <C 1 
AU. This is extremely improbable since the radius of the 
accretion disk should be even smaller, while typical sizes 
of accretion disks are more th an one order of magni tude 
larger. Is is worth noting that lAnglada et al.l (|2007l) ob- 
tain higher values of the proper motion velocities (~ 200 
km s _1 ), resulting in a smaller value for the estimated 
period, r p ~ 53 yr, that did not allow them to discard 
a precession-dominant scenario for the allowed range of 
masses. Thus, the wiggles observed in the plane of the 
sky have to be mainly caused by the orbital motion of 
the jet source. Precession, if any, is not important for 
shaping the jet and counterjet. 

In order to fit the shape of the jet and counterjet, we 
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Table 5 

Results obtained from the fit to the jet and counterjet knots. a 



Pnrfl tyi nfnr" 


(Units') 

1 \J 111 I 


Value 


Value 


Value 


Knots fitted 




ABCE b 


ABCDE b 


ABCDEZ b 


Initial rms 


(arcsec) 


0.360 


0.321 


0.272 


Fit rms 


(arcsec) 


n non 

u.uyu 


U.Uo / 


U. 1U4 


mill 


(M Q ) 


0.016 ± 0.006 


0.015 ±0.009 


0.014 ±0.006 


To 


(yr) 


110.1 ± 1.8 


110.6 ± 1.3 


113.9 ±1.7 


<Po C 


(rad) 


2.23 ±0.17 


2.15 ±0.17 


1.65 ±0.20 


To 


(AU) 


5.8 ±0.7 


5.7 ± 1.1 


5.7 ±0.9 


Vo 


(km s" 1 ) 


1.6 ±0.2 


1.5 ±0.3 


1.5 ± 0.2 


P d 


(deg) 











m j e 


(Mq) 


0.06-1.0 


0.06-1.0 


0.06-1.0 


Hc 




0.51-0.23 


0.50-0.23 


0.49-0.22 


m c { 


(Mq) 


0.06-0.30 


0.06-0.29 


0.06-0.29 


a 1 


(AU) 


11.4-25.1 


11.4-25.1 


11.5-25.6 



a Obtained for a distance D = 140 pc and assuming an inclination 
angle of the jet axis with respect to the plane of the sky i — 5°. 
b Knot velocities used: Vj = 98 km s _1 (A, B, C, E, Z4); Vj = 181 
km s" 1 (Zl, Z3, Z5); vj = 234 km s" 1 (D); Vj = 307 km s' 1 (Z2, 
Z6). 

c At the epoch t = 2010 Dec 01. 
d Adopted (see text). 

c Range of masses considered for the jet source (see text). 
f The first value is for the minimum value of rrij , and the second one 
is for the maximum. 

proceeded in various steps. For each step we used the 
procedure described in £15.61 First we fitted the "pure" 
orbital model to the knots of the jet that move at a simi- 
lar (low) velocity, and where the wiggling is more evident: 
knots A, B, C, and E. For these knots we adopted the 
same velocity, the error-weighted average of their veloc- 
ities, 98 km s^ 1 . The result of this first fit is shown in 
Table El 

The next step was to include in the fit the rest of knots 
of the jet, the high- velocity knots D. For these knots, 
we also adopted as velocity the error- weighted average of 
their velocities, 234 km s _1 (see Fig. [3]). The same values 
(within errors) of the parameters obtained for the first 
fit, were able to fit all the knots of the jet (see Table [5]). 

The third step was to include in the fit the knots of 
the counterjet. As stated before, the knots of the coun- 
terjet have very different velocities. We classified the 
counterjet knots into three groups. The medium- velocity 
group is formed by knots Zl, Z3, Z5a, and Z5b. Their 
error- weighted average velocity is 181 km s _1 . The low- 
velocity group contains only one knot, Z4, for which we 
adopted the same velocity as the low-velocity knots of 
the jet, 98 km s^ 1 . Finally, the extremely-high veloc- 
ity group, with knots Z2 and Z6, for which we adopted 
the error- weighted average of their velocities, 307 km s _1 
(see Fig. [3]) . Once again, the same values (within errors) 
of the parameters obtained for the first fit, were able to 
fit all the knots of the jet and counterjet (see Table [SJ. 

The last step was to include precession in the model, 
as a perturbation of the orbital motion (i.e. for /3 <C a). 
However, for this last fit the parameters of the precession 
were not well constrained by the data, and the rms fit 
residual did not improve with the inclusion of precession. 

Thus, we consider that the best fit to the shape of 
the jet and counterjet is an orbital model, with a mass 
function mul = 0.014±0.006 M Q , an orbital period r = 
114 ± 2 yr, and an orbital phase angle ip = 1.65 ± 0.20 



Table 6 

Favored set of parameters of the HH 30 binary system a . 



To* 

(yr) 


v c 
(km s- 1 ) 


(deg) 


c f 

rnj • 
(Mq) 


m c e,g 
(Mq) 


a e ' h 
(AU) 


114 ±2 


1.5 ±0.2 


95 ±11 


0.31 ± 0.04 


0.14 ±0.03 


18.0 ±0.6 



a Obtained for a distance D — 140 pc and assuming an incli- 
nation angle of the jet axis with respect to the plane of the 
sky i = 5°. 
b Orbital period. 
c Orbital velocity. 

d Orbital phase angle with respect to an axis perpendicular 
to the line of sight, at the epoch 2010 Dec 01. 
e Assuming a m ass of the system m = 0.45 ± 0.04 Mq 
ilPetv et al.li2"0M l. 
f Mass of the jet source. 
g Mass of the companion star. 
h Binary separation. 

rad at to — 2010 Dec 01. In Fig. [S] we show the rms 
fit residual as a function of the mass function m^ c and 
orbital period r , around the best-fit values. The shape 
of the jet obtained from this model is shown in Fig. [T0l 
where the continuum line indicating the shape of the jet 
has been obtained by interpolating the velocities adopted 
for the knots with a boxcar function, as shown in Fig. [3J 
From the fitted parameters we obtain a radius of the 
absolute orbit r = 5.7 ± 0.9 AU, and an orbital velocity 
v Q = 1.5 ± 0.2 km s . The orbital velocity obtained 
is low, similar to that obtained for other jets in young 
stell ar objects. For instance, in H H 111, v a = 3.1±1.9km 
s" 1 (jNoriega-Crespo et all [201 ll) and in HH 211, v D = 
1.6 ± 0.6 km s" 1 (|Lee et al.M20Toh . 

The values of the mass of the jet source and companion 
are not constrained by the fit. For the range of possible 
masses of the jet source rrij = 0.06-1 Mq (see above), we 
obtain that the companion has a range of masses m c = 
0.06-0.29 M , and the binary separation (radius of the 
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Figure 10. 2010 [SII] narrow-band image of the jet and counterjet of HH 30, between —30" and +60" of the HH 30 star. Left: Same 
scale in x and y. Right: x scale amplified a factor of 10. The measured positions of the knots arc marked with crosses. The rms x position 
of these knots with respect to the jet axis is 0'.'272. The continuum blue line is the jet shape calculated from the "pure" orbital model with 
a mass function mn% = 0.014 ± 0.006 Mq an orbital period r = 114 ± 2 yr, and an orbital phase angle ip = 1.65 ± 0.20 rad at to = 2010 
Dec 01. The rms fit residual in x is 0"119. The horizontal dashed lines mark the boundaries of the regions with different knot velocities. 
The knot velocities are indicated for each region. 



r elative orbit ) is in the range a = 11.5-25.6 AU. 

IPetv et ahl (|2006l ) find evidence of Keplerian rotation 
of the circumbinary disk in their CO observations with 
the PdBI, and derive a value of m — 0.45 ± 0.04 M Q for 
the central stellar mass. Using this value of m and the 
parameters of the fit, we obtain a mass rrtj = 0.31 ± 0.04 
M Q for the jet source, a mass m c = 0.14 ± 0.03 M Q for 
the companion, and a binary separation a — 18.0 ± 0.6 
AU. 

6. DISCUSSION 

Ang lada et al.l (|2007t ) carry out a detailed study of the 
morphology and proper motions of the HH 30 jet. They 
find that the proper motion velocities of the jet knots 
were of the order of 200 km s . They also find that 
the jet path shows a wiggling with a periodicity in the 
plane of the sky of 16" that translates into a period of 53 
yr using the derived proper motion velocity. It is shown 
that this period and the observed shape of the jet can be 



explained either by the orbital motion of the jet source 
or by precession of the jet axis because of tidal interac- 
tions of a companion star. In the first case, the binary 
separation would be 9-18 AU, while in the precession sce- 
nario the binary separation would be < 1 AU. In either of 
the two scenarios the binary separation is much smaller 
than the radius (~ 250 AU) of the disk observed w ith the 
HST (IBurrows et alfmM IStanelfeldt et alJl99flh . show- 
ing that it is a circumbinary rather than a circumstellar 
disk, and that it cannot play a role in the jet collimation. 
The true circumstellar disk associated with the HH 30 jet 
should be found at scales of a fraction of the binary sep- 
aration. In case of wiggling produced by orbital motion, 
the jet/counterjet system is expected to present plane 
("C") symmetry with respect to the equatorial plane, 
while in the case of precession the jet/counterjet system 
is expected to present point symmetry with respect to 
the positi on of the driving sour ce. Since the images ob- 
tained by lAnglada et al.l (|2007T > only cover a small por- 
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tion of the counterjet, they cannot discriminate between 
the two alternative scenarios. 

In this paper we presented new observations that al- 
lowed us to image and measure proper motions both 
in the jet and in the counterjet. In principle, these 
new data should be able to distinguish between the 
orbit al and precession scenario s. For instan ce, in HH 
111 ([Noriega-Crespo et"aTll201lD and HH 211 (jLee et al.l 
120101 ) plane symmetries in the jet/counterjet indicative 
of orbital motion of the driving source have been recently 
found. Unfortunately, the counterjet of HH 30 is weak, 
and appears less ordered than the jet, with larger varia- 
tions in the velocity of the knots. So, it is not expected 
to show a clear symmetry as a constant velocity wiggling 
jet would do. 

Interestingly, our new observations, covering a longer 
time span, allowed us to derive more accurate proper 
motions and reveal that the velocities of the knots 
are sm aller than the values estimated by lAnglada et all 
(2007). The new values of the proper motions imply that 
the wiggling period is about twice the previous value (114 
yr instead of 53 yr), allowing us to discard a significant 
contribution of precession to the observed wiggling of the 
jet (see $5/7]). 

Thus, based on our new observations we conclude that 
the observed shape of the jet is a consequence of the 
orbital motion of the jet source in a binary system. 
From our new results, and assuming that the central 
stellar mas s is m = 0.45 ± 0.04 M Q , as estimated by 
iPetv et ahl (|2006| ) from the CO kinematics of the disk, 
we obtain that the binary separation (i.e. the radius of 
the relative orbit ) shou ld be 18.0 ±0.6 AU. Interestingly, 
iGuilloteau et al.l (|2008l ) observe a hole of 37 ± 4 AU of 
radius in the CO and dust emission from the disk asso- 
ciated with HH 30. These authors interpret this hole as 
produced by tidal effects in a binary system, inferring a 
binary separation of 18 ± 2 AU, in excellent agreement 
with our result. Therefore, these results give strong sup- 
port to the binary interpretation, and in particular to 
the orbital motion scenario we proposed for HH 30. 

A binary separation of 18 AU corresponds to an an- 
gular separation of 0'.'13 at a distance of 140 pc. Thus, 
the proposed binary system associated with the HH 30 
jet can be angularly resolved by the new, expanded 
Jansky Very Large Array (JVLA) and by the Atacama 
Large Millimeter/submillimeter Array (ALMA). Unfor- 
tunately, according to our modeling of the orbital mo- 
tion, at the current epoch (~ 2012.5) the two stars are 
roughly aligned along the line of sight (0 O — 100°), re- 
sulting in a projected angular separation of 0'.'02, which 
is too small to be angularly resolved by the currently 
available instrumentation. We should wait four years 
until the projected angular separation will reach > 0'.'05 
and the presence of the proposed binary system could 
be confirmed by a direct observation of a double source 
and/or to measure orbital proper motions, as has been 
done previou sly in other binary you ng stellar objects (e.g. 
L1551-IRS5; iRodriguez et~alll200l 

When the HST observations revealed the silhouette of 
a flared edge-on disk with a radius of ~ 250 AU per- 
pendicular to the HH 30 jet, this system was considered 
the archetype of a jet/circumstellar accretion disk system 
in a young stellar object. However, the discovery that 
the central source is in fact a binary, and therefore, that 



the disk is circumbinary rather than circumstellar, has 
changed our understanding of this object. Tidal trun- 
cation of the disk in a binary system suggests that the 
radius of the circumstellar disk should be about ~ 1 /3 of 
the binary separation (e.g. iTerquem et"al~lll999h . Thus, 
we expect that the "true" circumstellar accretion disk 
associated with the driving source of the HH 30 jet must 
have a radius < 6 AU. Therefore, the search for this cir- 
cumstellar disk should be carried out at very small scales 
(~ 0'/05) 

Another interesting characteristic of the HH 30 sys- 
tem is that the jet and counterjet show clear differ- 
ences. The knots of the jet and of the counterjet can- 
not be grouped in pairs with similar separation from the 
star and similar velocity. HH 30 seems to be, in this 
aspect, completely different, for instance, from the jet 
and counterjet of HH 34, in which the jet and coun - 
terjet show a remarkable symmetry (Raga et al.l[20TlT ). 
Clear jet/counterjet asymmetries in the physical prop- 
erties have bee n found in several well studied jets, such 
as D G Tau B dPodio et alj|20ll. FS Tau B (jLiu et al.l 
[2012T). DC Tau (lAgra-Amboage et al.l l20H. or HH 30 
itself (jBacciotti et al.l I1999T ). Actually, strong velocity 
asymm etries between the red and blue lobes are fre- 
quent. iHirth et al.l (|1994j ) note that about half of bipo- 
lar jets present a strong asymmetry of about a factor 
two in velocity between the two lobes. In HH 30 we 
found that the average proper motion velocities in the 
counterjet are a factor of 1.7 higher than in the jet. We 
also found a higher velocity dispersion in the counterjet 
(see $3]). In the frame of stellar wind models, velocity 
asymmetries between the jet and counterjet have been 
attributed to differences in pressure. In the frame of 
magneto-hydrodynamical disk winds, jet asymmetries re- 
quire an asymmetry in magnetic lever arms or in launch 
radii between either sides of the disk, a situation that 
may occur naturally in an asymmetric ambient medium. 
For example, if the ambient radiation field is stronger 
on one side, this can increase the level of ionization on 
the surface of the disk on that side, leading to enhanced 
mass load on the magnetic field lines (smaller magnetic 
lever arm), and/or to a larger jet launching region (see 
iFerreira et all 12001 ) . 

However, HH 30 also shows differences in velocity in 
the same lobe of the outflow. For the jet, there are at 
least two velocities, a low velocity of ~ 100 km s _1 for 
knots A, B, C, and E, and a high velocity of ~ 240 km 
s _1 for knots D. For the counterjet we find knots close 
to each other with different velocities, and there are at 
least three different velocities, a low velocity of ~ 100 
km s _1 for knot Z4, a medium velocity of ~ 180 km 
s _1 for knots Zl, Z3, Z5a and Z5b, and an extremely- 
high velocity of > 300 km s _1 for knots Z2 and Z6 (see 
Sj3|. Indeed, we found that the knots Zl and Z2 (and 
perhaps also Z5 and Z6) of the counterjet appear to be 
launched nearly simultaneously with very different ve- 
locities. These strong differences in velocity for knots 
located in the same lobe of the outflow cannot be at- 
tributed only to asymmetries in the ambient medium 
between the two sides of the disk. Therefore, launch- 
ing from different radii in an extended disk wind sce- 
nario provides the simplest explan ation for the differ- 
ent ob served velocities. Following Blandfor d &: Pavnel 
(1982), the asymptotic jet velocity as a function of the 
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launch radius in the disk, n, is expected to be of the 
order of ~ 100(M*/0.5 Mq) 1 ' 2 ^/! AU)" 1 / 2 km s" 1 for 
a value of the magnetic lever arm parameter of ~ 10 
(i.e. "the extended warm disk wind scenario with mod- 
erate lever arms" de scribed, e.g. in iFerreira et al.l [20061 : 
iPanoglou et al.ll2"012L and references therein). Assuming 
that the jet source has a mass of 0.31 M Q (see above), 
we obtain that a range of launching radii n = 0.07-0.6 
AU is required to explain the velocity range of ~ 100- 
300 km s" 1 observed in the HH 30 jet/counterjet. This 
range of launching radii falls with in the range of values 
determined for T Tauri jets (e.g. [Anderson et al.1 120031 : 
iPesenti et al.1 120041 IFerreira et al.l 120061 and references 
therein), which range from 0.07 AU (the typical disk 
corotation radius) up to ~ 3 AU. 

7. CONCLUSIONS 

Using data from observations at six different epochs, 
we have measured the proper motions and studied the 
morphology of the HH 30 jet/counterjet system. Our 
main conclusions are summarized in the following. 

• The motion of most knots is essentially ballistic, 
but the jet/counterjet shows a large scale bending. 
This bending can be produced by a relative motion 
of the HH 30 star with respect to its surrounding 
environment of ~ 2 km s , caused either by a 
possible proper motion of the HH 30 star, or by 
the entraining of environment gas by the red lobe of 
the nearby L1551-IRS 5 outflow. Alternatively, the 
bending can be produced by the stellar wind from 
a nearby CTTS, identified in the 2MASS catalog 
as J04314418+181047. 

• The average velocity of the knots of the jet is about 
twice than that of the counterjet. Velocity differ- 
ences are also found in the same outflow lobe. Most 
of the knots of the jet move at velocities of ~ 100 
km s _1 , but several knots (knots D) move faster, 
at ~ 240 km s _1 . The velocities of the counterjet 
knots are less ordered, with knots moving, at least, 
at three different velocities of about 100, 180, and 
> 300 km s _1 . In particular, we identified at least 
a pair of knots in the counterjet (knots Zl and Z2) 
that apparently have been launched simultaneously 
with very different velocities. 

• The asymmetry in the average velocities of the jet 
and counterjet may be due to an asymmetric am- 
bient medium. The observed differences in velocity 
of the knots of the same outflow lobe can be at- 
tributed to different launching radii in an extended 
"warm" disk wind scenario. The range of launching 
radii in the circumstellar accretion disk required to 
account for the observed range of velocities of the 
knots is ~ 0.07 AU to ~ 0.6 AU, which falls within 
the range of typical values determined for T Tauri 
jets. 

• The observed wiggling of the HH 30 jet and coun- 
terjet is produced by the orbital motion of the jet 
source in a binary system. Precession of the accre- 
tion disk, if present at all, appears to be a minor 
contribution in shaping the jet. The best fit to the 



shape of the jet is obtained with a binary system 
with a mass function m/xj? = 0.014 ± 0.006 M©, an 
orbital period r D = 114 ± 2 yr, and an orbital ve- 
locity of v a = 1.5 ± 0.2 km s -1 . Adopting a value 
of the total stell ar mass of m — 0.45 ± 0.04 M & 
(jPetv et all 120061) . we obtain that the mass of the 
jet driving source is m = 0.31 ± 0.04 M©, the mass 
of the companion is m — 0.14 ± 0.03 M©, and the 
binary separation is 18 ± 0.6 AU. The same value 
of the binary separation is inferred from the size of 
the inner cavity observe d in the circumbinary disk 
(Guillote au et al.l 120081 ). assuming it is produced 
by tidal truncation of the disk in a binary system. 
This result strongly su pports our bina r y inte rpre- 
tation first proposed by Anglad a et al.l (2007). 



In summary, the HH 30 jet/counterjet system appears 
to originate likely from a binary system of two stars of 
masses 0.31 M© (the jet driving source) and 0.14 M©, 
separated 18 AU and orbiting with a period of 114 yr. 
The system is surrounded by a circumbinary disk with 
an outer radius of ~ 250 AU and with an inner radius 
~ 40 AU. The driving source of the jet/counterjet ap- 
pears to be surrounded by a circumstellar accretion disk 
with an outer radius < 6 AU, and probably with an inner 
radius of ~ 0.07 AU. The observed jet/counterjet knots 
reach velocities in the range ~ 100-300 km s^ 1 , arising 
from disk radii ranging from 0.07 AU to 0.6 AU. The two 
components of the binary system are currently roughly 
aligned along the line of sight, but their angular sepa- 
ration in the plane of the sky and orbital motions could 
be determined by future observations with the JVLA or 
ALMA 
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